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malization. It is consequently only necessary to calculate the running of 



D . the operators that do not vanish by the equation of motion. We show 



that our results are consistent with reparameterization invariance.Q 
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^ ; 1 Tree-level matching and class II operators 

Heavy quark effective field theory (HQEFT) is an approximation to QCD for quarks 
with masses that are large compared to the characteristic momentum scale of strong 
interactions. In the infinite- mass limit, it possesses Isgur-Wise spin-flavor symmetry, 
which facilitates the calculation of matrix elements of weak currents. The corrections 
to the theory for finite quark mass break Isgur-Wise symmetry, and these form an 
infinite series in powers of the reciprocal of the quark mass. 

There are two kinds of subleading operators which appear as terms in the cor- 
rected Lagrangian. Some operators, known as class I operators, do not vanish if 
the field is assumed to obey the leading-order classical equation of motion; others, 
known as class II operators, do vanish under this assumption. 



^ This is a long version of our paper, giving many details of the calculation. A somewhat 
shortened version (HUTP-96/A026) is being submitted for publication. 



Both types of operator appear in the Lagrangian obtained from matching to 
QCD. Tree-level matching with full QCD to order ^ yields the Lagrangian 

+^'^+A-^piD ■ vzp + {iD ■ (1) 



The leading-order equation of motion is 

D ■ vip+y = 0. 

Therefore, the operators 
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are class I, whereas operators such as 



Odv 
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(2) 



(3) 
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are class II. There is some freedom in defining what part of the Lagrangian is class 
I and what part is class II, since the class I operators may be defined to absorb part 
of the class II terms. Written in terms of (H), the Lagrangian is 



Ip+JD ■ V^p+^ + Okin + O 



mag 
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2m 



ij+^iiD ■ vfi;+^ + Oi + O2 



+^^+v[-m'^D -v-iD- v{ipf]i;+, + -^i^+A^D ■ vfi;+, (5) 
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In many treatments of HQEFT to order ^, Od v is simply thrown out of the La- 
grangian. It is not obvious, however, that simply throwing out the class II operators 
is correct at higher order in 1/m, and therefore we will discuss this point is some 
detail. 

The systematic method to remove redundant operators is via a suitable field 
redefinition. In the case (1), the field redefinition 
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Am 



32m2 



+v 



(6) 
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removes the class II operators; in terms of the redefined quark fields, it becomes 
-^0 + Okin + Omag + Oi + O2. Note that the field redefinition (Q) does not change 
the coefficients of the class I operators to order and so in fact is equivalent to 
the naive procedure of just dropping the class II terms. As we will show now, this 
is not the case in general. 

Class II operators can be removed from the effective Lagrangian to any desired 
order by use of the following iterative procedure. For this procedure to be possible, 
the the class I and class II operators must be defined such that they are separately 
Hermitian (note that this is not fulfilled by the operator definition in [Q). Let us 
assume that class II operators up to and including the order l/m"~^ have already 
been removed from the Lagrangian (n > 1, at the beginning of the iteration, no 
class II operators have been removed and n = 1). Collect the class II operators of 
order 1/m" in the form 

O class II = ;^^+-(^^ ■ + ■ ^)^+- 

i + i! 

Since Pj^ip+v = where P+ is the projection , this operator may be rewritten 

2 

O class II = ■ ^^+^ + ^^+^^ ■ ^)^+- 

Now perform the field redefinition 

^ (1 - (9) 

The P+ factor is there so that the identity ij'ip+v = may be satisfied by the 
redefined fields as well as the original fields. It is only necessary to include it if A 
is not written to commute with P+ to begin with. Clearly, this field redefinition 
removes class II operators from the Lagrangian up to and including 0(l/m"'). By 
iterating this process, class II operators can be removed to any desired order. 

Note that this field redefinition induces new terms at higher orders in the La- 
grangian. In general, we can split A^ in the form 

A"" = A1 + AliD ■ V (10) 

where A^ does not vanish when operating on an on-shell field with ivDip^^ = 0. If 
A" is non-zero, the field redefinition to remove the class II operators at order 1/m" 
will change the coefficients of the class I operators at order l/m"'*'^, and therefore 
in general it is not correct to simply throw out the class II operators from the 
Lagrangian. As A} = 0, this did not happen in the present case of but the 

naive method will yield incorrect results at order 
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2 Eliminating class II operators during renormal- 
ization 



In general, even after the class II operators have been eliminated by field redefini- 
tion, a renormalization of the theory will induce class II terms. The theory with 
class I operators alone cannot be renormalized without a revision of what is meant 
by a renormalization. Therefore, when calculating without class II operators, an 
infinitesimal renormalization must be followed by an infinitesimal field redefinition 
that removes the class II operators. Then, in general, the Lagrangian will possess 
the form 



L = ilpj^^D ■ V1p+^ + CkinOkin + CmagOraag + CiOi + C2O2, (U) 



to order \ 



The renormalization of the order ^ operators at one loop is well known [^, 



Ckin is not renormalized at all, and 



d_ _ _%2 



B - n—r - ^ r 

yraag " A ^mag / . \2^mag 



(4vr) 



^C^ag hiSUiN). (12) 



The coefficient of O^.y is also renormalized at one loop. However, the renormal- 
ization is entirely multiplicative, so if this operator is eliminated from the Lagrangian 
by a field redefinition after tree- level matching, there will be no Od v term after an 
infinitesimal renormalization either. 

C.L.Y. Lee attempted to renormalize the order ^ operators, but Lee's di- 
vision of operators into class I and class II parts was not Hermitian, so it was not 
actually possible to perform the field redefinition necessary to eliminate the class 
II terms. Here we renormalize the Lagrangian to order ^ and one loop using the 
Hermitian class I operator definitions (^. Also, unlike Lee, we define coefficients 
only for the local operators in the Lagrangian and derive j3 functions in terms of 
those coefficients. In 0, S. Balk et al. discuss the renormalization of a variant 
of HQET constructed from a sequence of Foldy-Wouthuysen transformations to 
0{l/w?). While the operator Oi is not present in their basis, their results seem to 
be consistent with ours regarding the running of O2. For a comparision with other, 
more recent calculations |^ see section |^. 

Since the class II terms induced by renormalization are all of order -K, when 
renormalizing the Lagrangian, there is no need to actually calculate the renormal- 
ization of all of the class II terms and calculate the field redefinition necessary to 
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remove them during renormalization. The effect of the field redefinition on the La- 
grangian to order ^ will be equivalent to simply throwing out the class II terms 
induced by renormalization. 

The class I operators induce various Feynman vertices, shown in Figures [l| and 
0. The Feynman rules are given in the appendix. Under background field gauge, as 
described in Abbott explicit gauge invariance establishes relations between the 
coefficients of the various vertex terms that are preserved under renormalization; the 
vertices arising from a single gauge- invariant operator are multiplied by the same 
running coupling constant. 

3 Renormalizing O2 to one loop 

Because of this explicit gauge invariance, in order to renormalize an operator, it is 
only necessary to examine the part of each IPI loop diagram's divergent part which 
has the form of the simplest of the operator's Feynman vertices. Furthermore, if the 
vertex factor in question has multiple terms, it is only necessary to look for one of 
them, provided that it is not possible to produce the sought-after term using class 
II operators. This immensely simplifies the task of renormalizing the operators. 

O2 is the easier of the order operators to renormalize. 02's one-gluon ver- 
tex has a term which is rather difficult to create in diagrams involving the other 
operators: 

^k'^p^a-^v'^T^. (13) 

This is a natural thing to look for, because it puts strict constraints on the form 
of a divergence, so few diagrams will contribute to its renormalization. For any 
divergence to renormalize this interaction, it must not only have the same tensor 
structure, but also depend on the external momenta k' and p in the correct way. 
They must both be contracted with a"^, rather than with v, with themselves, or 
with each other. 

3.1 Diagrams with one internal gluon line 

It is worthwhile to consider a very common situation in which these criteria are 
usually not satisfied. Suppose that a one-loop diagram has only one internal gluon 
line. It is perfectly legitimate to label the momenta so that the loop integration 
variable, q, is the momentum of that internal gluon. Then, no factors of p and no 
factors of k' will appear in the gluon propagator. 

There will, in general, be terms with factors of k' and p in the denominators of the 
heavy quark propagators. However, there a momentum always appears contracted 
with the heavy quark's four- velocity v. Therefore, expanding these propagators in 
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powers of p will never yield a term in which k' or p is contracted with anything other 
than V. 

When doing a loop integral, the typical procedure is to combine denominators 
using some sort of Feynman parameter, then shift the integration variable so as to 
make the Euclidean integral hyperspherically symmetric. In this case, though, all 
the propagators but one are heavy quark propagators. The shift in the integra- 
tion variable only transforms q into q — Xv, where A is a Feynman parameter with 
dimensions of mass. Factors of q do not turn into factors of k' or p later in the 
calculation. 

This would seem to indicate that in a diagram with only one internal gluon line, 
any factor of (|13D in the divergence must explicitly appear in the product of the 



vertex factors before the loop integral is done. However, there is another possible 
complication 0. Integrating over all momentum space can transform products of 
loop momenta into factors of the metric. In particular, if /(g^) is any scalar function 
that depends on q^ only via g^, 

(14) 

When considering which diagrams have the correct contractions to contribute to 
the renormalization of O2, it is necessary to take into account situations in which 
things that ought to be contracted with each other are both contracted with the 
loop momentum. 

Fortunately, the terms we are looking for already have two factors of the external 
momenta in them, so this mechanism requires the presence of at least four factors 
of momentum in the numerator. Therefore, it only operates in a few diagrams. (In 
fact, it turns out not to lead to any new diagrams with one internal gluon line in 
the renormalization of O2, because there is no way to get four factors of momentum 
from the requisite sets of vertices.) 

These considerations allow the elimination of many diagrams with little effort. 



3.2 Diagrams with two internal gluon lines 

The vast majority of relevant diagrams therefore contain two internal gluon lines, and 
a three-gluon QCD vertex connected to the external gluon leg. In these diagrams, 
factors of p may arise from places other than the various vertices involving quarks, p 
appears in the three-gluon QCD vertex. It also appears in at least one of the internal 
gluon momenta, which makes it easier to obtain it in the quark-gluon vertex factors. 

It is still not necessary to worry about factors of p in the propagator denominators 
when renormalizing O2, since at worst they will produce factors of p ■ v and p^, not 
factors of p contracted with the a tensor. 
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Any factors of k' must still appear explicitly in the quark-gluon vertex factors 
prior to loop integration. It is possible to label the momenta in the loop so that k' 
only goes through heavy quark lines. In denominators it is contracted only with f , 
and it never appears in the expression for the shift in q after combining denominators. 

However, now the identity (|l^ is important, because it is possible to obtain four 
factors of momentum with the use of the three- gluon QCD vertex in addition to 
the heavy quark vertices. It is necessary to consider both diagrams in which (13) 
appears explicitly in the numerator prior to loop integration, and diagrams in which 
some of the factors in ([13|) are contracted with q^q^ rather than with each other. 

These facts allow the elimination of some of the diagrams with two internal 
gluons. Diagrams with more internal gluons do not appear, because they have 
either too many loops or too many legs. 



3.3 Operator insertions that do not contribute 

The term (|T3|) is spin-dependent, so it is not necessary to consider diagrams con- 
taining only spin-independent vertices. Double insertions of Okini and insertions of 
Oi, will therefore not contribute to the running of O2. Either Omag or O2 itself has 
to be in the diagram somewhere. 

Double insertions of Omag cannot produce (|1^) either, because Omag has no 
factors of k' in any of its vertices. Therefore, the only contributions to the renor- 
malization of O2 must come from diagrams involving an insertion of Okin o,nd Omag, 
or diagrams with a single insertion of O2. 

In this case, there is also no need to subtract out divergences of the form of the 
class II operators, because none of the class II operators has a one-gluon vertex with 
a term of the same form as (|T3p . The only spin-dependent class II operator of order 

■\ is -ip,^{iD ■ v,G"^aaH}ip+v, and the heavy quark residual momentum k' is 

never contracted with a in the vertices. 



3.4 Diagrams with Okin and Omag 

There are many one-loop diagrams involving one Okin vertex and one Omag vertex; 
fortunately, most of them do not contribute to the term in question. 

As said above, in all cases the correct factor of k' must arise explicitly from the 
vertex factors, and they must be contracted as in (0), except that factors of q^q'^ 
may appear in place of the metric. 

If the one-gluon vertex of Okin appears, that vertex may not be connected to the 
external gluon leg or to a leading-order quark- gluon vertex, because in neither case 
will k' end up contracted with the a tensor in Omag, or with the loop momentum q. 
It must be connected to an internal gluon line. 
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If the gluonless Okin vertex appears on an internal quark line carrying some linear 
combination of k', p, and q, the only way (1T3|) may be obtained is in the term with 
a factor of k' ■ q, so that (|I^) can leave k' contracted with a. 

In either case, the k' term in the Okin vertex has no factor of p. Therefore, a 
factor of p must appear somewhere else. If there is only one internal gluon line, then 
the factor of p must come from the one-gluon Omag vertex, by the general argument 
above. But it cannot do so. If the momentum in the Omag vertex is p, then that 
vertex must be connected to the external leg, and there is no way to contract the 
factor of a with k' or q. 

The diagram must, therefore, have two internal gluon lines. If k' or q is to end up 
contracted with the a tensor, an internal gluon line has to contract with Omag- There 
are only four such diagrams, shown in Figure ^. The calculation of the contribution 
of Figures ^-b is described in more detail than the other parts of the calculation; 
it is a typical example of the process of calculating these divergences. 

The amplitudes corresponding to these diagrams are 



d'^q 



gfabc[gp^i~'2pp) + gp^{2q - p)^ + g^,p{2p^ 



{p — qYq'^{k' — q) ■ V 



from the diagram in Figure and 



I 



X 



gfabc[gpni-2pp) - ^p«(2g - p)p + g^p{2p^) 
—i 



(p — qYq^{k' — p + q) ■ V 



X 



from Figure The projection operator in the heavy quark propagator may be 
ignored, since all operators are assumed to be sandwiched between heavy quark 
fields, and in both cases the projector commutes with one of the vertex factors. 

Taking only the terms with a single factor of k' contracted with a, and adding 
the two diagrams together, gives 

_ 2m^^Tak' a"" [ A qa{^q-p)p f _±_ ^ i \ (^7) 

{2my ° ^ J {2tt)'^ {p — qYq^ \<i-v {q — p) ■ v J 

where N is the number of colors. The quadratic factors in the denominator may be 
combined using the usual method of Feynman parameters; shifting the integral to 
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make it symmetric, and extracting only the term linear in p yields 

-2/V.p^T..>„.-/^^ (18) 

The heavy quark propagator factors remaining in the denominator may now be 
combined with the rest using the usual (for HQEFT) identity 

1 _ (n + m-1)! 2™A™-MA 
. v)"" ~ {n - l)!(m - 1)! Jo {q^ + 2Xq ■ ^ ' 

where A is a sort of Feynman parameter with dimensions of mass. Then the integral 
over momentum space is 



oo 



Again shift the integral, q ^ q — \v, 

(2^^g2_p)3 (21) 

and remove the term that is odd in g; now the integral over A is easy, and the 
amplitude becomes 



Under DRMS in 4 — e dimensions, the amplitude is 

(47r)2e 4m2^"^^^ ^^^^ 

In Figures ^-d, contributions of the form of (|T^ arise from factors of the form 
k' ■ qqaPi3<y°'^ 1 via ([T^). Taken together, the relevant part of their contribution is 

After using Feynman parameters to combine the denominators and making the 
integrand symmetric, ( p!^ may be used to obtain a term of the desired form, which 
is 

^^■i^>.-"VT.. (25) 
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3.5 Diagrams with O2 

Some of the diagrams containing an insertion of O2 may be eliminated in an analo- 
gous manner to the ones considered above. Here it is not necessary to take (|14D into 
consideration, because an O2 vertex and a three-gluon QCD vertex can together 
contribute at most three factors of momentum in the numerator. 

The three-gluon O2 vertex is completely nonderivative, so there is no way to 
obtain the desired factor of k' from that vertex. 

The term in the two-gluon vertex that has the k' lacks any other factors of 
momentum, so the factor of p must come from somewhere else. Therefore, any such 
diagram must have two internal gluon lines. The only possibility is Figure |^a. 

With the one-gluon vertex, there will be no contribution from diagrams in which 
O2 is connected to the sole internal gluon line, since then the gluon momentum in 
the vertex is g, and there are no other vertices with momenta in them. However, 
there is also the option of connecting O2 directly to the external gluon leg, so that 
a factor of the desired form comes directly from the O2 vertex. That is Figure 

There are also diagrams with two internal gluon lines, one of which contracts 
with O2S one-gluon vertex; these are Figure ^ and Figure ^d. 

The calculations proceed much as before. In each diagram, the terms to look for 
are the ones in which the momenta are appropriately contracted, keeping in mind 
that factors of a loop momentum q can turn into factors of p when the integration 
variable is shifted, if p flows through an internal gluon line (but not if it only flows 
through an internal heavy quark line). The important term contributed by Figure 
|a is 

The i is the loop integral's symmetry factor. Combining denominators in the usual 
way reveals that the p in the denominator contributes only terms of higher than 
linear order in p and may therefore be neglected (the story would have been different, 
had there been terms with factors of q in the numerator). The divergent term linear 
in p is 

-i^-A^>-"VT. (27) 



Figure ^ makes a small contribution, because of a group-theoretic factor of -2^; 
the term with the correct factors of k and p is 

^ Sl'.)^>.-"V/^-5^, (28) 



(2m)2 V 2N V ""^^ J {27T)^q^{q-v 
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which comes to 

Figures ^ and ^ contribute in much the manner of Figures |^-b. The term of the 
desired form in Figure ^jc, after combining the gluon denominators, shifting the loop 
variable, and doing the integral over the Feynman parameter, is 

- f / ^^^^ (30) 

Figure ^'s contribution is identical. Adding the two together, combining with the 
heavy quark denominators and performing the requisite integrals gives 

+ (47r)2e 4m2^"^^^ ^'^^ ^^^^ 
3.6 The (5 function of 02's coefficient 

These, then, are all of the divergent pieces of T021 the part of the IPI three-point 
function with the same form as (18). Adding them all together, including the coeffi- 
cients from the Lagrangian, and putting in the contribution from the single vertex. 



02 



Co 



9 
Aim? 



2 \ „ 1 



{An)' 



In^ 



k'^Ppa'^'^v^'Ta + ■■■ (32) 



where the dots represent convergent terms not dependent on /i. Now the j3 function 
for C2 may be determined by solving the renormalization group equation for Tq^'- 

[3kin is zero {Ckin does not run and is just equal to 1), the [3mag term is higher order 
in g than the others, and explicit gauge invariance in background field gauge means 
that the terms with I3g and 7^ cancel. The anomalous dimension of the heavy quark 
field is 

Solving for to order (7^, 

P2{g{fi), Ckin, CmagifJ')) = J^;^ {4:N C kinC mag) (35) 
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Remarkably, there is no multiplicative renormalization of C2', the diagrams with O2 
insertions are completely canceled by the term from the heavy quark anomalous 
dimension. At the scale where matching occurs, Cmag = Ckin = C2 = I, and 

p2ig{mQ), Ckrn, Cmag{mQ)) = (36) 

Since every diagram except ^ contains a three-gluon QCD coupling, it is also 

interesting to consider the case of a U{1) gauge theory, in which is the only IPI 

1 

diagram that renormalizes O2. Then ^ lacks the group-theoretic factor of that 

is present in the SU{N) case, and the heavy quark anomalous dimension is j^^- 
The solution of the RGE to order is 

/^2 = ^(+4-4)C, = (37) 

so, for f/(l), C2 does not run at all at one loop. 



4 Renormalizing Oi to one loop 
4.1 Eliminating class II terms 

The renormalization of Oi is slightly more involved, conceptually and mathemat- 
ically, for two reasons. First, there are more diagrams to consider; second, the 
simplest vertex of Oi has no term that cannot be produced by class II operators 
as well, so it is necessary to calculate the part of each diagram that renormalizes 
one of the class II operators, and subtract it out. Fortunately, these two difficulties 
cancel each other out to some extent, because many of the extra diagrams turn out 
to renormalize only the class II operator. 

The term in the one-gluon vertex of Oi that gives the least trouble is 

(pV)T,. (38) 



This term also appears in the one-gluon vertex of the class II operator -ip, {iD ■ 

V, D'^}ip^^. However, in this operator the factor p'^v^ always shows up in the com- 
bination {—2k' ■ p + p'^)v'^, and the k' ■ pv'^ term is not produced by any other local 
operator to order Therefore, to subtract out the renormalization of the class II 
operator, all that is necessary is to add, in the divergent term of each diagram, half 
the coefficient of 

^.{k'.pv^Ta (39) 
8m^ 
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to the coefficient of (|38|). 

Both of these terms are spin-independent. Therefore, contributions to the renor- 
mahzation may come from double insertions of Okin, from double insertions of Omag, 
or from single insertions of Oi itself. 



4.2 Diagrams with one internal gluon line 

The expression ( pHD has no factor of k' in it, so it is now possible to obtain nonzero 
contributions to the running of Oi from diagrams in which no factors of k' appear in 
the vertices. There may be a contribution if the numerator, prior to loop integration, 
contains factors like k' ■ p, p^, k' ■ qp ■ q, or p ■ qp ■ q. If there is only one internal 
gluon line, these factors cannot arise from the propagators, for precisely the same 
reasons as in the O2 case. Routing the loop momentum q through the lone gluon 
line reveals that its denominator is simply g^, and factors of k', p, or q in the quark 
denominators are contracted with v. 

Therefore, in such diagrams, the factors of k' and p in the divergence must appear 
explicitly in the vertex factors prior to loop integration, and the external momenta 
must be contracted with each other or with q'^q'^. 

Class II contributions may also be neglected. If the factors of the form of and 
k' ■ p appear in the combination — 2k' ■ p, then the divergence just renormalizes a 
class II operator, and it need not be considered. 



4.3 Diagrams with two internal gluon lines 

When there are two internal gluon lines, additional complications arise. Factors 
of can and do arise in the gluon propagators. Therefore, when calculating the 



divergent term of the form of (|38|) , it is necessary to expand the loop integrand 
in powers of and retain the zeroth- and ffist-order terms. For this reason, the 
calculation of these diagrams is more involved than for Oi. 

4.4 Diagrams with Okin 

Terms like (^) and (|39D arise from many diagrams with two Okin vertices. 



4.4.1 The gluonless vertex 

In particular, it seems as if contributions to both ought to come from any diagram 
with a no-gluon vertex on a quark line carrying some linear combination of k' and p 
(and possibly the loop momentum q) . This is indeed the case for the diagrams with 
two internal gluon lines (as in Figures ^o-d). 



13 



If there is only one internal gluon line, the possibilities are more restricted. If 
there is also only one gluonless vertex, then there can be at most three factors of 
momentum and cannot provide the desired contractions. Routing q through 
the gluon line implies that the quark line on which the Okin vertex sits may only 
carry either fc'±g or k' —p±q. In the former case, the gluonless Okin vertex provides 
neither a term like ( pSD nor one like (|39|). In the latter case, the terms appear in the 
class II combination and may be ignored. 

If there is one internal gluon line but two gluonless Okm vertices, it becomes 
possible to obtain contributions via (|14|) . If both of them are on the same side of 
the external gluon leg, however, the momentum factors are of the same form as 
those of a single vertex, and the diagrams may be excluded by the arguments in 
the previous paragraph. The only class I contribution occurs when the gluonless 
vertices are on opposite sides of the external gluon leg, as in Figure 



4.4.2 The one-gluon vertex 

The one-gluon vertex has uncontracted momenta in it; factors of k' ■ p and p'^ may 
be produced by putting it in a diagram with two internal gluon lines (Figures ||a-c) 
or by contracting it with another one-gluon Okm vertex (Figure ^). 



4.4.3 Eliminating the two-gluon vertex 

The two-gluon vertex is completely nonderivative, so any diagram with that one in 
it would have to get its factors oik' -p or p^ from somewhere else. If there is only one 
internal gluon line, then, as stated above, these factors would have to come from the 
other Okin vertex. They can't, since they could only arise from the no-gluon vertex, 
and that only supplies them in a class II combination if there is only one internal 
gluon line. Therefore, the two-gluon Okin vertex can only contribute if there are two 
internal gluon lines. 

But if the diagram is to be IPI, then one of those has to connect to something 
other than the two-gluon Okin vertex, which leaves it with an extra leg. So this 
vertex doesn't contribute at all. 

The diagrams in Figure |^ are therefore all of the contributions to the renormal- 
ization of Oi from double insertions of Okin- When calculating them, it is necessary 
to do the expansions in p^ mentioned above. 

For example, consider Figure §a. After the usual manipulations to combine the 
gluon denominators, shift the integrand, and throw out terms proportional to /c'^ or 
k' ■ V OI to four- vectors other than or g^, the leftover amplitude looks like 



(2m)2 2 7o J i2TiY[q^ + x{l-x)p'^Y{q-v) 
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X [{x^ + x+l)p^ + + {l + 2x)p-q-4k' ■q-2{l + 2x)k' -p] (40) 
Expanding the factor with the denominators 

1 1 2x(l-x)p2 



[g2+x(l-x)p2]2 (g2)2 (^2)3 



+ ••• (41) 



reveals the presence of additional terms in the result. The relevant divergent 
terms turn out to be 



-.3 



, ( ANk' ■ p + / , , L (42) 

(2m)2 ''V 2 ^ / y (27r)4 (g2)2 V ; 

Subtracting out the class II contribution, by adding the coefficient of the p^ term to 
half the coefficient of the k' ■ p term, makes this 

i5^^r„pV + (//) (43) 

where (II) refers to the part that renormalizes an operator that vanishes by the 
equation of motion. 

Figures and |^ together contribute 

- (^^"(*' /IfiW (^^) 

which is 

f^^T-P^.." + (//). (45) 

and Figure ^ contributes 

p^T.(-8iV.'.p + 2iVp^).„./|^^ (46) 

which is 

(47r)2e Sm^''^^'^ ^^''>- ^^'> 
Figures ^-f have a different group-theoretic structure; |^'s contribution is 

16 (yf^/i" ig 
~ iV (47r)2e8m^ 

which is 



rT.fc ■ pv^" (4J 



^'^^ '^T,pV + (//). (49) 



A^(47r)2e8m2' 
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and fig:kinkinonef's is 



1 \ ^ [ d'^q q^{p^ -2k' -p) -Ak' -qp-q 
-Ta / TTTZU :JTZ-Xi (50) 



2N J 4m2 J (27r)4 q^{q ■ v)'^ 

which, after combining denominators, shifting the integration variable, and applying 



(|T^), becomes 

8 g^n' ig 



3N (47r)2e8m2 



r,pV + (JJ). (51) 



In this case the entire class I contribution comes from (]1^) . 



4.5 Diagrams with Omag 

Double insertions of Omag will produce no factors of k' ■ p, since in such diagrams 
the only factors of k' are in heavy quark propagators. Therefore, it is only necessary 
to seek out factors of p'^v'^. 

If the diagram has only one internal gluon line, at least one of the Omag vertices 
will be connected to it, so its vertex will contain a factor of q rather than p. Then 
there is no way of obtaining two factors of p explicitly in the quark-gluon vertex 
factors; ([T^) does not apply, since there is no way to obtain four factors of momentum 
in the numerator. 

There must be two internal gluon lines. That leaves only Figure |^. 

In this calculation, the projector in the heavy quark propagator must not be 
neglected, since it does not commute with a'^'^ 0. When extracting the spin- 
independent part of this diagram's amplitude, one may use the identities 

^af, I l±i j ^ - 2vf'v^ + spin-dependent terms (52) 



and 



'l+f 



<ypvp'^p'^ = -gf.up'^ + 9tiv{p-vf + p^p^ + p^Vf^Vy 

—v^PuP ■ V — PfjiVuP ■ V + spin-dependent terms (53) 



where both expressions are assumed to be sandwiched between heavy-quark spinors, 
so that factors of ^ may be dropped at the beginning or end of a term. 
Figure ^s contribution is 

^^T„p^... (54) 

3(47r)^e 8m^ 

There is no class II term. 
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4.6 Diagrams with Oi 

In the multiplicative renormalization from insertions of Oi, no k' -pv^ terms arise, so 
there are no class II renormalizations to subtract away. Again, ( p^lf ) does not apply 
since there are at most three factors of momentum in the numerator. 

If there is only one internal gluon line and it contracts with any of the Oi 
vertices, a factor of cannot arise explicitly in the vertex factors, since the only 
vertex with two factors of momentum in it is the one-gluon vertex, and in that case 
the momentum will be q. The only diagram left that has one internal gluon line is 
Figure 0b. 

The other possibilities have two internal gluon lines, one of which contracts with 
an Oi vertex. These are Figures |^, |^c, and When calculating these diagrams, 
it is necessary to expand combined gluon denominators as series in p^, as before. 

Figure ^a's symmetry factor is "2"; cLS for Figure Ha. The relevant divergent term 

is 



(47r)2e 9>m? 

Figure 0b gives 



Tap'v^'. (55) 



much in the manner of Figure ^d's contribution to the renormalization of O2. Finally, 
Figures 0c and 0d together give 

Ta—V / TTT^ (57) 
our ~ 

or 



8m2 " 6 ^ J (27r 



3(47r)2e 8m' 



4.7 The function of Oi's coefficient 

The calculation of (5\ is analogous to that of (i^- The three-point IPI function of the 
form of (pHD, with class II divergences subtracted, is 



To, = {C^ 



X -^pVr, + --- (59) 
8m^ 

where the dots represent convergent terms not dependent on Solving the RGE 



~ 3n) + -^^rna, + + ^ ) ^^1 



^ In/i 



(An) 
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to order gives 

/3l(5'(/^),C'fcin,Ci(/i)) 



9' 



(47r) 



-6A^ + 



16 

3iV 



^2 



fern 



lOiV ,4iV 

2 '^mag + ~3~^ 



C'fcin = 1. At the scale where matching occurs, Ci 



a 



mag 



1, and 



N+ — ] 



3Nj (47r)2' 



(61) 



(62) 



Figures ^-f and ^ are the only diagrams present in an abelian gauge theory. 
1^ contributes to the renormalization of Oi in exactly the same way that Figure ^ 
renormalizes O2. |^-f's contribution is as above only without the SU{N) factor of 
Therefore, at one loop in a U{1) theory. 



Pi 



9' 



(47r)2 
32 

y 



32 

-yC^L + (4-4)Ci 
9' 



(63) 



Other than the gauge coupling constant and the field normalizations, this is the only 
thing that runs at one loop to order -\ in a U{1) theory. 



5 Reparameterization invariance is satisfied 

In HQEFT, the division of the quark momentum into a large part mv^ and a small 
residual momentum is arbitrary, as long as the residual momentum remains 
small. Therefore, there exists a symmetry of HQEFT called reparameterization 
invariance, in which the four-velocity and residual momentum change so as to leave 
the combination mv^ + invariant. ||10|| This places constraints on the coefficients 
of the correction terms. 

However, it is first necessary to define how the heavy quark field transforms 
under a reparameterization. There are at least two different forms of reparameteri- 
zation invariance in the literature, that of Luke and Manohar and that of Chen 
TT[|. Straightforward calculations show that the Lagrangian obtained from tree- 



level matching (before the field redefinition that removes the class II operators) is 
invariant under Chen's transformation, but not under Luke and Manohar's. Chen's 
transformation 



2 



1 + 



2m + iD ■ V 



(64) 
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does the following to the Lagrangian £ = ip^^Aip^^ to order ^: 

6C = (^6A + [A, ^m6vx] + {^^ f } + + As/-^^ + O (-1) 

(65) 

where pj_ = p — D ■ v^. (Because of the term proportional to m, we would need to 
know about ^ terms in the Lagrangian to evaluate what the transformation does 
to order -K.) 

Applying (^5|) to (|TT]) induces a change in the Lagrangian 

bL = '?/' . ^(1 - Ckin)iD ■ 5v4)+v + - — -ip^^ [-D ■ 5vD ■ v - D ■ vD ■ 5v 

4m 

+ (1 - 2C^ag + C2){iD^^5v''a^,D ■ v - iD ■ vD^dv^'a^,)] (66) 

The field-redefinition procedure we have used is equivalent in its effect on the 
running of the class I operators to throwing out all class II terms. Insertions of the 
class II terms does not affect the running of the class I terms, because the class I 
parts of the loop diagrams involving class II operators are finite. The poles that arise 
from solutions of the classical equation of motion are eliminated by the vanishing of 
the class II vertices. Therefore, our procedure will yield the same running for the 
class I operators that we would get if we kept all of the operators, class I and class 
II. 

The tree-level Lagrangian, including the class II operators, is symmetric under 
the RPI symmetry in (|65|). Therefore, even using our procedure, as long as the 
regularization procedure preserves the RPI symmetry in (^), the renormalized La- 
grangian to order must be symmetric under this transformation, up to terms 
resulting from the action of the transformation on the removed class II operators. 
The —D ■ 6vD ■ v — D ■ vD ■ 6v term in ( |66D results from the action of the transfor- 



mation on Oo v, so there is no reason to expect that term to vanish. The others, 
however, cannot be so obtained. Therefore, reparameterization invariance sets the 
constraints 

2C^ag = C2 + I (67) 

At matching, all of these constants are equal to 1, and the constraints are satisfied. 
Under running, in order to maintain these relations, it must be that 

2/3mag = /?2 (68) 

It is well known that the first relation holds to the orders that have been studied. 
The second also holds at one loop, according to our results. The running satisfies 
reparameterization invariance in the form described by Chen. 
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6 Comparison with other recent calculations 



In this section, we compare to results in 0, |^ 

While this work was in preparation, a paper by Balzereit and Ohl was posted on 
the net |Q, which also calculates the renormalization of the l/w? operator. Their 
technique is quite different from ours; most importantly, they retain all class II 
operators, including Od.^. However, insertions of class II operators do not induce 
any class I counterterms, since the amplitudes resulting from class II insertions are 
finite. The poles that arise from solutions of the classical equations of motion are 
eliminated by the vanishing of the class II vertices. Therefore, a calculation using 
our technique, which assumes that the Lagrangian contains only class I operators, 
will yield the same j3 functions (or, equivalently, anomalous dimensions of local 
operators and time-ordered products) obtained in for the class I operators. 

Balzereit and Ohl use a slightly different operator basis, but their basis for the 
class I operators is the same as ours up to class II terms and a sign difference in Oi. 

Even given the class II terms by which our basis for the order \ local class I 
operators differs from [Q, the results still have to be equivalent. With both operator 
bases, our technique can be used and the calculation will differ in no essential respect. 
For example, to renormalize Oi in the basis of [Q, one could look for terms of the 
form 2k' ■ pv'^, and subtract out class II terms of the form {2k' ■ p — p'^)v'^ by adding 
twice the coefficient of 2p to the coefficient of 2k' ■ pv^. This is manifestly the 
same calculation, and so the results must be the same, because the class II terms by 
which our operator bases differ do not affect the calculated coefficients of the p'^v'^ 
and k' ■ pv^ terms. 

Indeed, our results are equivalent to those of [Q. They define Oi with the 
opposite sign, which reverses the signs of the contributions to Oi from double Okm 
and double Omag insertions. Our definitions of (3 functions also differ from their 
definitions of anomalous dimensions by a further factor of —'^g^/ (47r)^. Taking these 
into account, the terms in our (3 functions are equivalent to the various elements in 
their anomalous dimension matrix ^^pflyg- The nonzero anomalous dimensions in the 
first column correspond to the coefficients of Ci, C^j„, and C^^^ in our expression 
for j3i. The nonzero anomalous dimension in the second column corresponds to the 
coefficient of CkinCmag in our (32- 

Very shortly before we posted this paper to the net, another calculation of the 
renormalization of the order operators by B. Blok et al. appeared. This 

calculation, like that of Balzereit and Ohl, retaines all class II operators, and ex- 
presses results in the form of anomalous dimensions mixing local operators with 
time-ordered products. The class I operator basis in is the same as ours up to 
total derivatives, except for a sign difference in the definition of Omag- The results 
for the class I operators also agree with ours, when this sign difference is taken into 
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account; the first two columns of the matrix (22a) in |Q agree with the terms in our 
beta functions for the case N = 3. 



7 Conclusions 

We have calculated the running of the heavy quark effective field theory Lagrangian 
to order using a technique in which continuous field redefinition removes opera- 
tors from the Lagrangian which vanish according to the classical equation of motion. 
Our results are consistent with symmetry under the reparameterization transforma- 
tion of Our results are inconsistent with those of and agree with those of 
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A HQEFT Feynman rules to order ^ 

The operators of HQEFT induce various Feynman vertices. In addition to the rules 
hsted below, there are of course the usual QCD rules for gluons; we use background 
field gauge with a Feynman-like gauge prescription (at one loop, the gauge fixing 
parameter does not run and Feynman-like gauge is OK), so the gluon rules are as 
given in Abbott ||^ with a = 1. The covariant derivative and gluon field-strength 
tensor are defined as follows: 



It is convenient to express the rules in terms of the outgoing heavy quark's residual 
momentum k''^ (taken as flowing out of the vertex) and the gluon momenta p'^, q'^, 
and r'^ (which flow into the vertex), a, b, and c are external gluon color indices. The 
vertices from the subleading operators are pictured in Figures |I] and |^. 

A.l Leading-order Lagrangian 
A. 1.1 propagator 



(69) 



(70) 




(71) 



A. 1.2 



one-gluon vertex 



igTaV 



(72) 



A.2 



kin 



A.2.1 



no-gluon vertex 



/2 



(73) 



2m 
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one-gluon vertex 



two-gluon vertex 



Omag 

one-gluon vertex 



two-gluon vertex 



^^9 uy f rp 

-p. — Cr JabdJ-d 

2m 



one-gluon vertex 



two-gluon vertex 



8m^ 



[(? - P) ■ vg'"" - (p + 2qrv'' + {2p + qYv^lUaTd 



8m? 

three-gluon vertex 



• 3 

[{g^'v" - g^^v'')ncdfdaeTe + {g'^'V' - g'^^V^fcadfdbeTe 



8m^ 

+ {g'^'v'^ - g'-'^vnfabdfdceT, 

02 

one-gluon vertex 

^ {2k'^Pfsa"f'v''-p-v{2k'-p)^a^^)Ta 



8m2 
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A. 5. 2 two-gluon vertex 



8m2' 

+{-p^a''^v^-q^a''''v'' + {p - q) ■ vaniT^T^^}] (82) 
A. 5. 3 three-gluon vertex 

+ {a''^v^ + a^'v'')hM{Tc,Ta}] (83) 
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k' 
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